Shape dynamics during deposit of simple metal clusters on rare gas matrices 
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Using a combined quantum mechanical/classical method, we study the collisions of small Na 
clusters on large Ar clusters as a model for cluster deposit. We work out basic mechanisms by 
systematic variation of collision energy, system sizes, and orientations. The soft Ar material is 
found to serve as an extremely efficient shock absorber. The coUisional energy is quickly transfered 
at first impact and the Na clusters are always captured by the Ar surface. The distribution of the 
collision energy into the Ar system proceeds very fast with velocity of sound. The relaxation of 
shapes goes at a slower pace using times of several ps. It produces a substantial rearrangement of 
the Ar system while the Na cluster remains rather robust. 
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I. INTRODUCTION 

Clusters on surfaces have been much investigated dur- 
ing the past decade as can be seen, e.g., from the sequence 
of recent ISSPIC proceedingaii^i^i^. The topic remains of 
great interest, especially in relation to the synthesis of 
nano-structured surfaces. A possible way is here the di- 
rect deposition of size selected clusters on a substrate^i^. 
When deposited on a surface, a cluster undergoes a sig- 
nificant modification of its electronic structure and ionic 
geometry, because of the interface energy, the electronic 
band structure of the substrate, and the surface corruga- 
tion. These questions have been widely investigated from 
the structural point of view both at the experimental- 
and theoreticaliSiiiii^iiliiiii^ii&iii sides. The theoretical 
description of deposition dynamics is nevertheless a very 
demanding task. A fully microscopic approach treating 
all atoms and electrons in detail is extremely involved and 
applicable only to very small samples if at all. Most theo- 
retical descriptions thus employ simple molecular dynam- 
ics approaches with effective atom-atom forces. There 
are, however, situations where a detailed description of 
electronic degrees of freedom is desirable, for example 
when metal clusters are involved where electronic shell 
effects are known to play a role in forming the structure. 

A step forward is to combine a fully microscopical 
treatment of the cluster with a much simplified descrip- 
tion of the surface. This is a valid and efficient com- 
promise for inert substrates as, e.g., for the deposit of a 
metal cluster on an insulator surface. Such an approach 
was explored, e.g., in the case of Na clusters on NaCl 
jj^jis Density Functional Theory was used for 

the Na electrons coupled to the surface via an effective in- 
terface potential, itself tuned to ab-initio calculations^^. 
This approach, however, ignores any excitation and/or 
rearrangement of the surface itself, which supposes ex- 
tremely inert materials. The next step is obviously to 
restore a minimum of surface degrees of freedom, still 
having electronically inert substrate in mind. Such a 
model was recently developed for the case of Na clusters 
embedded in rare gasesiSi^iiS where the cluster is again 



treated in full detail and the rare-gas atoms through clas- 
sical dynamics of position and dipole polarization. Such 
a "hierarchical" approach is justified for electronically in- 
ert substrates. It has much in common with the coupled 
quantum-mechanical with molecular-mechanical method 
(QM/MM) often used in bio-chemistry22i2^i2i. The now 
explicit, although simple, handling of substrate atoms 
still restricts calculations to finite systems. The model 
developed i n^o, 21^25 ^ nevertheless, allows to explore a suf- 
ficiently large range of sizes to see the appearance of 
generic behaviors on the way towards the bulk. 

The goal of this paper is to study the dynamics of de- 
position of a small Na cluster on a finite model of an 
Ar surface. To that end, we consider Nag, Nay, or Nas 
clusters as projectile on Ar^r clusters of various sizes {N 
between 7 and 87) as target. We analyze the kinetics of 
deposition in terms of the energy transfer between cluster 
and substrate, a point for which the explicit dynamical 
surface degrees of freedom are crucial. We shall also an- 
alyze the way in which the metal cluster adapts to the 
substrate and how the substrate gives way. This is done 
by tracking the evolution of the shapes of both partners 
during the deposition process and by varying Ar cluster 
size as well as initial kinetic energy. For this purpose 
we shall compare several small Na clusters with differ- 
ent shapes in the initial state. It should finally be noted 
that the considered setup (Na cluster on Arjv) has two 
aspects. First, it is an example of cluster-cluster colli- 
sions involving partners of very different size and nature, 
as such a topic of its own interest. Second, we have in 
mind an exploratory study of deposition on a Ar surface 
for which we consider large, but finite systems, a method 
which has been used also in experimental studies, see 
e.gj^. The widely used alternative is to simulate the in- 
finite surface in terms of a small, but periodically copied, 
simulation box, see e.g;^^. This better allows to include 
long range effects but limits the structural rearrangement 
of the surface. Thus modeling the surface by finite clus- 
ters is probably more appropriate for our intention to 
analyze the shape dynamics. 

The paper is organized as follows. Section [Til provides 
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a short presentation of the model and of the systems 
considered in the test cases. In section IIIIl we analyze 
the deposition dynamics in terms of trends and energies. 
In section IIVI we discuss the shape dynamics in more 
detail. 



II. MODEL 

As stated above we use a hierarchical approach in 
which the metal cluster is treated in full microscopic de- 
tail while the "surface" is described at a classical level 
granting each atom mobility as a whole and a dynamical 
dipole response. This is justified by the clear hierarchy 
of electronic binding in metals vs. rare gases. The many 
different ingredients make a complete description of the 
model rather bulky. We give here a short account of the 
approach and refer to^i^ for a detailed layout. 

The Na cluster is treated using quantum-mechanical 
single-particle wavefunctions {(p„(r, t), n = 1, . . . , N^}, 
for the valence electrons coupled non-adiabatically to 
classical molecular dynamics (MD) for the positions of 
the Na ions {R/,/ = l,...,7VNa}. The electronic 
wavefunctions are propagated by time-dependent local- 
density approximation (TDLDA). The electron- ion in- 
teraction in the cluster is described by soft, local pseudo- 
potentials. This TDLDA-MD has been validated for lin- 
ear and non linear dynamics of free metal clusters^^'^^. 

Two classical degrees of freedom are associated with 
each Ar atom: center-of-mass {Ra, a = 1, . . . , iVAr}, and 
electrical dipole moment {da, a = 1, . . . , iVAr}. With the 
atomic dipoles, we explicitely treat the dynamical polar- 
izability of the atoms through polarization potentials^. 
Smooth, Gaussian charge distributions are used for Ar 
ionic cores and electron clouds in order to regularize the 
singular dipole interaction. The Na+ ions of the metal 
cluster interact with the Ar dipoles predominantly by 
their charge. The small dipole polarizability of the Na"*" 
core is neglected. The cluster electrons do also couple 
naturally to the Coulomb field generated by the atoms. 

The polarization potentials describe the long-range 
Coulomb part of the interactions. There remains to ac- 
count for the short-range repulsion. The repulsive Na^- 
Ar potential is taken from^i. The pseudo-potential for 
the electron- Ar core repulsion has been modeled in the 
form of'^^, slightly readjusted by a final fine-tuning to the 
NaAr molecule (bond length, binding energy, and opti- 
cal excitation). For the atom- atom interactions, we use 
a standard Lennard- Jones potential. A Van der Waals 
interaction is added, computed via the variance of dipole 
operatorsi2i2ii2^. It provides a contribution to the long- 
range interaction which is crucial to produce the faint 
binding of the ground states for small Na-Ar systems. 
Its dynamical effects stay merely at a quantitative level. 

We use standard methods^ for the numerical so- 
lutions. We solve the (time-dependent) Kohn-Sham 
equations for the cluster electrons on a grid in coor- 
dinate space, using time-splitting for the time propa- 



gation and accelerated gradient iterations for the sta- 
tionary solution. In the present calculations, we fur- 
thermore use the cylindrically-averaged pseudo-potential 
scheme (CAPS) as an approximation for the electronic 
wavefunctions'^'^''^^. Both the Na ions as well as Ar cores 
and valence clouds are treated by classical molecular dy- 
namics (MD) in full 3D. We have checked that a full 3D 
treatment of the electronic wavefunctions leads to simi- 
lar results for the chosen systems which stay all close to 
axial symmetry. 



System 


Na-Na dimer 


bulk Ar-Ar 


Na-Ar dimer 


Bond energy 


800 meV 


50 meV 


5 meV 



TABLE I: Bond energies of dimers entering the calculations 
discussed in the text. 



Table U shows the binding energies for the three pos- 
sible dimer combinations. There is a clear hierarchy of 
binding where the Na cluster is strongest bound, the Ar 
cluster one order of magnitude less bound, and the Na-Ar 
binding another order of magnitude smaller. This leads 
us to expect a fragile attachment of the two clusters to 
each other while the clusters as such remain intact with 
some readjustment of shape, particularly at the Ar site. 




FIG. 1: Initial configurations for the deposition of Nae (white 
balls) on Ary (left), Ar43 (middle), and Argy (right). 

Three Ar clusters of different sizes are used as "sub- 
strate". At J, Ar43, and Ar87. The structures of these 
clusters are optimized for the given model by simulated 
annealing. The collision partners are various Na clusters 
taken in their ground state structure, again obtained by 
simulated annealing but this time coupled to static Kohn- 
Sham iterations for the electrons. Most of the results pre- 
sented below have been obtained for the strongly oblate 
Nae cluster, composed of a ring of five ions and an outer 
ion, as can be seen from Fig. [TJ Top configuration refers 
to the case where the pentagon hits first the Ar cluster. 
Reverse configuration corresponds to a Nag up-side down. 
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We also consider a few other clusters of comparable size 
but different shapes, especially Nay and Nag. In Nay, 
the central pentagon is complemented by one top and, 
symmetrically, one bottom ion, which altogether delivers 
a much less oblate shape than Nag. The Nag cluster, 
in turn, is built as two layers of four ions. Mind it is 
electronically close to sphericity because it contains the 
magic number of eight electrons. And indeed, the shape 
of free Nag is, as far as it can be for a small finite cluster, 
close to spherical. In all cases, we take care of orienting 
the Arjv relative to Nag in such a way that it presents 
to the Na cluster the largest possible planar "surface". 
This can be seen in Fig. [T] where the initial configuration 
for the collision between Nag and Ary, Ar43 and Argy are 
displayed : The metal cluster initially faces respectively 
1, 5, and 4 Ar atoms. Note that in the largest system, the 
top Na ion starts above an interstitial position, whereas 
it is initially placed above an Ar atom in the two other 
cases. 

The dynamics is started by giving the system a rela- 
tive boost, with a "substrate" at K temperature. We 
have studied the effect of " substrate" temperature up to 
50 K. No significant differences are observed what the 
deposition mechanisms and associated energy transfers 
are concerned, although the highest temperature repre- 
sents a thermal energy of the same order of magnitude 
(about a factor 2-3 smaller) than the kinetic energy in 
Na. This shows that relevant effects are primarily of 
potential nature, especially the balance between short 
range (Ar core) repulsion and long range (Ar polarizabil- 
ity) attraction, both aspects well taken into account in 
our model. Note also that experimentally, deposition of 
metal clusters on rare gas substrates are performed at 
very low temperature (typically below 25 K)^^^. Us- 
ing just K temperature is thus fully sufficient for the 
present aim of a qualitative study. 



III. ANALYSIS OF DEPOSITION MECHANISM 
A. A visual example 

We first consider in Fig. [2] a typical deposition scenario 
for the example of Nag on Ar43 with initial kinetic en- 
ergy £^kinO = 13.6 meV per Na atom. The Nag cluster is 
accelerated when coming closer to the Ar cluster, a long 
range polarization effect which is counterbalanced by the 
core repulsion only at very short distance. The pinning 
process proceeds stepwise with several bounces before the 
metal cluster is finally attached to the rare gas cluster. 
The overall shape of the Na cluster remains basically un- 
changed while the Ar cluster is somewhat rearranged, as 
was expected from the relative binding energies in Ta- 
ble |T1 It should be noted that the initial kinetic energy 
per atom of the impinging cluster is smaller than the 
bonding energies of both Na2 and Ar2 dimers but about 
2.5 larger than the binding energy of the NaAr dimer. 
Not surprisingly neither pure Na nor pure Ar bonds are 



broken, but in spite of the faintness of the NaAr bond 
one observes asymptotic stitching. This implies that an 
efficient kinetic energy transfer has taken place to cool 
down the Na-Ar interface, the excess kinetic energy be- 
ing transferred to internal energy (potential and heat) of 
the Ar and/or Na clusters. A part of the potential en- 
ergy is visibly used to reshape the Ar surface, in order 
to provide a more convenient contact plane for the Nag 
cluster (see late times in the figure). The energy trans- 
fer mechanisms may a priori depend both on the initial 
kinetic energy and on the Ar cluster size. This will be 
analyzed in the coming sections. 





0.5 ps 


1.5 ps 


2.5 ps 


4.25 ps 


7.25 ps 
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^^^^^^^^^^^^^^^^^^ 


^^^^^^^^^^^^ 









FIG. 2: Snapshots of the deposition of Nas on a planar site 
of Ar43, with an initial kinetic energy of 13.6 meV/ion. Time 
slots are indicated on each panel. 



B. Systematics on Ar cluster size 

We now turn to the question of the influence of the 
Ar cluster size and its heat capacity on the deposition 
process. Fig. [3] shows results for Nag deposited on vari- 
ous Ar clusters all for the same initial kinetic energy of 
-Ekino = 68 meV per Na atom. The separation between 
the centers of mass of the two cluster is initially 30 Qq 
for all systems. The left panels show the detailed z coor- 
dinates (symmetry axis and direction of collision) for Na 
ions and Ar atoms. It is obvious that, at this projectile 
energy, the Ary cluster is broken into pieces after the im- 
pact while both Ar43 and Argy are massive enough to re- 
sist the impact. Thereby the smaller Ar43 is visibly more 
perturbed (molten?) while the heaviest sample Argy even 
maintains its shell structure. The kinetic energies shown 
in the right panels explain that trend. Almost the same 
total kinetic energy of about 0.3 eV is transferred in all 
three cases, independent from the Arjv size. Distribut- 
ing that equally over the Ar atoms associates a typical 
temperature to the Ar cluster of about 200 K for Ary, 
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30 K for Ar43, and 10 K for Arg7, which fits perfectly to 
the observation of break-up, melting and stability. The 
largest sample thus serves as a reasonable model for de- 
posit on a surface in this exploratory study. However true 
bulk is, of course, thermally very inert due to its large 
calorific capacity and one would expect even smaller val- 
ues of temperature in the case of bulk deposition under 
the same kinematic conditions. 
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coordinates of the Ar atoms in the left panels below 2 ps. 
The perturbation propagates like a straight line through 
the Ar layers with a speed of 20-30 ao/ps. Long time 
scales, on the other hand, remain for the final relaxation 
processes. The successfully captured Nae cluster contin- 
ues to bounce and oscillate with a relaxation rate of order 
of 10 ps. The fine-tuning of the Ar cluster shape seems 
to proceed at about similar slow scale. Note, however, 
that this shape evolution initiated by the sound wave 
might differ for an infinite substrate where the sound 
wave would dissolve into deep layers. In the small finite 
systems here, the wave is dissipated by diffuse scattering 
from the opposite surfaces of the Ar cluster which turns 
collective energy into heat. A large part of the Ar kinetic 
energy seen in Fig. [3] thus becomes thermal energy. The 
case of largest Ar cluster, Argy, behave somewhat differ- 
ent. It presents cleaner surfaces as compared to smaller 
clusters and it shows a faint reverberation of the wave 
when the latter reaches the edge of the cluster, as is ob- 
served in the lower left panel of Fig. [31 The reflected 
wave bounces back with almost the same velocity to the 
side which qas hit initially. There it returns some mo- 
mentum to the Nag, as can be seen in the small revival of 
its kinetic energy at about 3 ps (lower right panel). This 
reverbation is, of course, due to the finite size of ArgT 
and would not appear in the case of an infinite substrate. 
Nevertheless, from the metal cluster point of view, this 
reflection seems to change very little its final shape and 
distance to the Ar "surface" . In that sense, this system 
still mimics the gross features of deposition on an infinite 
surface. 



FIG. 3: Collision of Nae (thick lines) with an initial kinetic 
energy of 68 meV/ion, on Arjv (thin curves) for N = 7 (top), 
= 43 (middle) and iV = 87 (bottom). Left panels: z- 
coordinates; right panels: total kinetic energies, as a function 
of time. 



A closer look to the kinetic energies in Fig. [3] shows a 
few interesting differences. The common feature is an ini- 
tial energy transfer which leaves about half of the initial 
energy, i.e. around 0.3 eV, as kinetic energy in the Ar and 
very little kinetic energy for the Na cluster. The other 
half is invested in potential energy to provide the large 
spatial rearrangements. Differences are seen in further 
evolution. For Ary all rearrangement is finished after ini- 
tial encounter and the kinetic energies remain unchanged 
later on. The severe perturbation of the medium sized 
Ar43 leads to ongoing rearrangements with long lasting 
exchange between potential and kinetic energy of the Ar 
cluster. As a side effect, there occurs also some re-heating 
of the Nae cluster. Smaller reflow of potential into kinetic 
energy is seen for the heaviest Argr. 

Fig.[3]also spans very different time scales. The energy 
transfer at the moment of impact is surprisingly fast. The 
Na kinetic energy is fully given up within less than 0.5 
ps and the energy distribution within the Ar cluster pro- 
ceeds as a sound wave. This can be seen from the z 



C. Systematics on initial kinetic energy 

It is expected that the deposition depends on the velo- 
city of the projectile. We analyze this effect in the case of 
Nag deposition on Argy, with varying E'kino- Results are 
plotted in Fig.[l]for i^kino between 6.8 meV and 136 meV 
per Na atom. The striking feature is the overall similar- 
ity of all the different cases. This concerns in particular 
the very fast and almost complete energy transfer from 
the Na kinetic energy to the other degrees of freedom. A 
common feature is also the quick capture of the Na clus- 
ter and its long standing residual oscillations. Of course 
details of the scenarios differ from one case to the next. 
The perturbation of the Ar cluster increases naturally 
with increasing initial energy. There is also an interest- 
ing effect at initial times. The kinetic energy of the Na 
cluster first increases before contact. That is due to the 
medium range polarization interaction which is attrac- 
tive. The additional acceleration depends on the initial 
velocity. More energy is gained for the slower velocities 
because the cluster moves for a longer time in the attrac- 
tive regime. Thus the kinetic energy is almost quadrupled 
for the lowest i?kinO = 6.8 meV while only 10% effect is 
left over for the fastest collision here. Besides these sub- 
tle differences in detail, the energy deposit proceeds in all 
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cases the same way, namely in an extremely quick fashion 
leaving basically no residual kinetic energy at the side of 
the Na cluster projectile. One thus finds that in a rel- 
atively large range of projectile kinetic energies, the Na 
projectile is glued to the Ar target. Except again for 
oscillations in the most energetic case, the overall struc- 
ture and position (with respect to the Ar target) of the 
deposited Nae seems to depend very little on the initial 
energy. 
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FIG. 4; ^-coordinates (left) and total kinetic energies (right), 
as a function of time, in the collision of Nae (thick curves) on 
ArgT (thin curves), for four different initial kinetic energies 

-EkinO- 



IV. TIME EVOLUTION OF SHAPES 

A. Cluster distance 

The above calculations show that the relaxation times 
of the whole deposition process goes beyond times we 
can reasonably simulate. However, we have seen from 
Figs. [3] and [4] that around 6-7 ps the kinetic energy of 
the Ar cluster and the cluster distance seems to reach a 
constant mean value. The closeness of the two clusters 
can be quantified in terms of a "deposit coordinate" d 
that we define as 



A^NaA^Ar |R/ " Ra 



-1/4 



(1) 



where the index / (respectively a) refers to the Na ion 
(respectively Ar atom) cores and R/ (respectively Rq) 
to their position. The inverse puts weight on the closest 
partners. The actual power 4 in Eq.(IT|) is a matter of 
decision: we calculated with powers n = 1 up to n = 



Eui„n=136 meV 



7 and found that n = 4 provides a good compromise 
between averaging and details. 
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FIG. 5: Deposit coordinate d, defined in Eq.([T}, as a func- 
tion of time, for Na6@Ar43 (dots and dashes) and NagQArgr 
(full lines), with -EkinO = 13.6 meV/ion (left panel) and 136 
meV/ion (right panel). Results in both systems are presented 
for the top (dashes and thin curve) and the reverse (dots and 
thick curve) configuration of the Nag. 



Fig. [5] shows the coordinate d for two different sizes 
of the Ar clusters and for two different i?kinOi i-^. initial 
Na kinetic energies per ion. The pattern are surprisingly 
similar for both energies and also similar for both sizes. 
The Na cluster is catched at first impact. There remain 
some oscillations in distance of about 2-3 ag which relax 
very slowly (order 10 ps). There is a global difference to 
the extend that the distance is closer for Na6@Ar43. This 
is due to the larger rearrangements in that case which, 
in turn, provide a better contact area for the Nag. 

We have also varied the orientation of the Na cluster. 
The standard configuration for Nae was such that the 
ring of five Na ions stood closer to the surface while the 
top ion was pointing away (the situation denoted as "top" 
in the figure). Within axial symmetry, there is also the 
reverse situation where the top ion is facing towards the 
surface (denoted "rev" in the figure). The result for re- 
verse initial configurations are also shown in Fig. (5] They 
are almost indistinguishable from the top configuration. 
This indicates that the Nag remains basically unchanged 
in all cases. This will be corroborated in the next section 
studying shape dynamics. 



B. Global moments 

The results discussed up to now indicate that there 
may be substantial rearrangements at the side of the Ar 
substrate. Less is yet known about the shape of the Na 
cluster: do we have a soft landing, some plastic deforma- 
tion or a wetting behavior? In order to quantify these 
questions we have performed a shape analysis in terms of 
the first three multipole moments both for Nag and Ar^r. 
These moments are given by 
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where p is either the number of Na atoms A^Na or the 
number of Ar atoms A^Ar, and x, y and z are the coordi- 
nates of the Na (Ar) atom with respect to the center of 
mass of the Na (Ar) cluster. The r.m.s. radius r stands 
for the overall extension (monopole moment) and the de- 
formations are parameterized as dimensionless quantities 
which have immediate geometrical meaning independent 
of system size. For example, a value of \[32\ ~ 0.8 is a 
large quadrupole deformation with axis ratio of about 
2:1. 
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FIG. 6: Three first multipole moments (r^), P2 and Ps, as 
a function of time, for Nae (thin dashes and lines) deposited 
on Ar43 (thick dots) and Arg? (thick full lines) for -EkinO = 68 
meV/ion. 



Fig. [5] shows the three momenta for the Na and Ar 
subsystem in the cases Na6@Ar43 and NagQArgy for the 
moderate initial kinetic energy -Ekino = 68 meV per 
Na ion. The shape of Nag is rather rigid in any case. 
There are some deformation oscillations short after im- 
pact which relax within about 3 ps. These oscillations 
are predominantly caused by the outer ion. The ring 
is tightly bound and stays more robust. Note that the 
relaxation is much faster than for the overall bouncing 
oscillations of the cluster (see Fig. [S|). That is related 
to the binding properties shown in Table HI that is, the 



NaAr binding is softer. The Ar clusters, after the impact 
with Nag, increase slightly in size due to the heating by 
energy absorption. The growth is relatively larger for 
the smaller Ar43 which acquires a higher temperature 
as discussed above. The Ar clusters do also undergo a 
strong persistent change in deformation towards a size- 
able oblate and somewhat pear-like shape. They obvi- 
ously accommodate their configuration as to establish a 
most compact combined system. 

For the case of Nag@Ar43, the effect of initial kinetic 
energy and orientation (Na cluster in "top" or "rev" con- 
figuration, see section HV A[) are shown in Fig. [T] For the 
large Skino = 680 meV/ion (right panels), the Ar43 emits 
rapidly several atoms, thus yielding a rapidly increasing 
radius and octupole moment while a more moderate evo- 
lution emerges for the less violent i^kino = 136 meV/ion. 
The collision with the Nag top configuration is more vi- 
olent because the whole ring bounces first and at once. 
This gives rise to an even higher Ar deformation for both 
initial energies. However, at the side of Nag, mean value 
and amplitude of the shape oscillations are similar in all 
cases, for both configurations and for both i?kinO- Note 
the change of sign of for the reverse configuration. 
This indicates that the outer ion oscillates through the 
pentagon, within the given time once for i?kinO = 136 
meV and twice at 680 meV. But the other two moments 
stay very robust in spite of the violence unload upon the 
Ar43 subsystem. 



C. Wetting behavior; comparison with Nay and 

Nag 

We thus find that the Nag is very robust under any con- 
ditions when deposited on the rather soft Ar material, ft 
maintains its oblate shape and overall radius. The top 
ion is more loosely bound and may undergo larger oscil- 
lations which become apparent in the octupole momen- 
tum. The results for Nag thus suggest that wetting of 
the surface is rather unlikely in that combination of ma- 
terials. However, Nag is not a very conclusive test case 
as it has already a close to planar structure. One ion still 
sticking out of the surface plane looks not so dramatic. 
In order to countercheck, we considered two neighbor- 
ing clusters, namely Na7 (less oblate than Nag) and the 
almost spherical Nag. The results for the multipole mo- 
ments are shown in Fig.[51 They are very similar to what 
we have seen already for collisions with Nag. The Nay 
or Nag clusters maintain basically their shape while the 
Ar basis undergoes a substantial rearrangement to cre- 
ate space for landing. There is not the slightest hint of a 
wetting. There are, of course, small shape oscillations in 
the Na clusters due to the internal excitation at the time 
of impact. It is interesting to note that these oscillations 
differ in detail between the two clusters. The Nag with 
its magic electron number 8 seems to be more rigid and 
shows generally less oscillation amplitude. 
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FIG. 7: Three first multipolc moments, as a function of time, 
for Nae (thin curves) deposited on Ar43 (thick curves), with 
-Bkino = 136 meV/ion (left) and 680 meV/ion (right). On 
each panel are compared the moments obtained from the top 
(dots and dashes) and reverse (full lines) configurations. 
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CONCLUSION 



In this paper wc have investigated slow collisions of a 
small Na cluster with a large Ar cluster as a test case for 
deposit on an Ar surface. To that end, we have employed 
a hierarchical model treating the Na cluster in full detail 
by time-dependent density-functional theory and the Ar 
atoms at a thoroughly classical level with position and 
polarization as dynamical variables. We have studied in 
detail the influence of various parameters on the deposi- 
tion scenario, namely initial collision energy, Ar cluster 
size, and relative orientation of the Na cluster. 

We have found that, except for the case of the smallest 
Ar system, the basic scenario is robust in the sense that 
it depends very little on the collision parameters. The Na 
clusters are well bound while Ar binding is much softer 
and the Na-Ar binding even less. As a consequence, the 
Ar substrate constitutes a very efficient shock absorber. 
The impinging Na cluster is stopped immediately at first 
impact and attached to the Ar system while the latter 
takes over practically all initial Na energy. It seems im- 
possible to tunc conditions under which the Na cluster 
is reflected. Enhancing the initial energy rather leads to 
destruction of the Ar system. The initial large energy 
transfer from the Na cluster to the Ar system is very 
quick, taking less than 0.5 ps. The transferred energy 
is distributed also very quickly over all Ar atoms, prop- 
agating like a shock wave with speed of sound through 
the medium. About half of the energy goes into potential 
energy and part of it is used up for a large spatial rear- 
rangement of the Ar cluster which aims to provide a most 
compact combined system. These rearrangements pro- 
ceed on a slower time scale of 4-10 ps. At the side of the 
Na cluster, we have observed that in all cases its shape is 
little affected by the dynamical processes during deposit. 
It always safely attaches to the surface but only loosely 
bound and keeps nearly its free shape. The interface in- 
teraction is too wc^ak to induce a wetting mechanism for 
this material combination. 

Calculations with larger clusters and planar surfaces 
for the same combination confirm the above findings. It 
is now most interesting to check other combinations of 
materials. Work in that direction is in progress. 
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